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Final  Report;  Algorithms  for  Point  Set  Congruence 
ARO  Grant  No.  DAAL03-92-G-0378 

Under  the  support  of  this  grant,  the  investigator  studied  the  question  of 
congruence  between  two  sets  of  points.  The  problem  draws  its  motivation 
from  computer  vision.  One  point  set  can  be  thought  of  as  a  model,  and 
the  other  as  an  image.  A  given  model  could  represent  a  symbol,  such  as  a 
letter  of  the  alphabet,  or  a  landscape  or  other  environment.  An  image  is  that 
which  is  seen  by  an  observer.  We  wish  to  know  if  an  observed  image  matches 
a  certain  model. 

Before  the  tenure  of  this  grant  began,  this  investigator  had  begun  research 
into  several  formulations  of  this  problem.  The  investigator  had  improved 
upon  results  of  other  researchers  by  introducing  a  technique  that  perturbs 
the  input  points  slightly  in  order  to  exploit  the  resulting  geometric  structure 
of  the  point  sets.  The  investigator  had  applied  this  technique  to  the  case 
of  equ2d  cardinality,  planar  point  sets.  The  resulting  algorithms  were  called 
“approximate  algorithms,”  because  they  traded  some  precision  for  faster  run¬ 
times. 

The  primary  purpose  of  the  project  supported  by  this  grant  wais  to  ex¬ 
pand  research  in  point  matching,  by  applying  the  perturbation  technique  and 
the  approximate  algorithm  paradigm  to  more  complex  problem  formulations. 
During  the  tenure  of  the  grant  the  investigator  has  successfully  worked  to¬ 
wards  this  goal.  The  main  results  of  this  research  are  contained  in  the  paper 
“Generalized  approximate  algorithms  for  point  set  congruence”.  This  paper 
is  still  undergoing  revisions,  and  a  copy  of  the  most  recent  draft  is  included 
in  this  report.  The  paper  introduces  a  device  called  the  (e,  fc)-map,  which 
is  a  more  general  measure  of  point  set  congruence  than  the  one  previously 
examined  by  the  investigator.  The  paper  discusses  ways  to  construct  the 
(e,  fc)-map.  While  the  investigator’s  previous  work  was  limited  to  equal  car¬ 
dinality,  planar  point  sets,  this  current  work  has  expanded  attention  both 
to  point  sets  of  unequal  cardinality  and  to  point  sets  in  higher  dimensions. 
The  investigator  has  even  begun  the  study  of  projective  congruence,  in  which 
one  compares  a  2-dimensional  image  with  the  family  of  urojections  of  a  3- 
dimensional  model. 
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1  Equal  cardinality,  planar  point  sets 

The  case  of  equal  cardinality,  planar  point  sets  was  the  main  focus  of  this 
investigator  before  the  tenure  of  this  grant.  The  major  results  of  this  work 
are  contained  in  the  paper  “Approximate  decision  algorithms  for  point  set 
congruence,”  written  by  this  investigator  with  Stefan  Schirra.  Final  revisions 
on  the  paper  were  made  by  the  authors  during  the  period  of  the  grant,  and 
the  paper  has  been  accepted  in  revised  form  by  the  journal  Computational 
Geometry;  Theory  and  Applications.  The  content  of  this  paper  is  important 
to  the  main  body  of  work  pursued  under  the  grant,  since  the  recent  work 
has  studied  generalizations  of  the  problem  formulated  in  this  paper.  We 
summarize  below  the  problem  studied  in  “Approximate  decision  algorithms 
for  point  set  congruence,”  and  the  approximate  algorithm  paradigm. 

The  paper  studies  the  problem  of  e- congruence.  Assume  that  we  are  given 
two  planar  point  sets  A  and  B  of  equal  cardinality,  n,  and  a  tolerance  value 
e.  An  isometry  is  a  mapping  in  R.^  that  preserves  distances.  We  say  that  A 
and  B  are  e- congruent  if  there  exists  an  isometry  /  and  a  bijection  £  :  A  B 
such  that  dist{t{a),I(a))  <  e,  for  all  a  €  A,  where  dist{‘,-)  is  the  distance 
function  of  our  chosen  metric  (typically  £,2)-  Intuitively,  we  can  approach 
this  problem  by  thinking  of  each  point  set  being  put  on  a  piece  of  paper; 
the  sets  are  e-congruent  if  it  is  possible  to  slide  one  piece  of  paper  over  the 
other  and  match  the  points  in  such  a  way  that  the  points  of  each  pair  are 
within  distance  e  of  each  other.  Important  specieil  cases  of  this  formulation 
are  obtained  by  restricting  the  isometry  to  be  a  translation,  T,  or  a  rotation 
around  a  fixed  point  d,  Rd. 

Determining  if  point  sets  A  and  B  are  c-congruent  for  a  given  e  is  referred 
to  as  the  decision  problem.  The  optimization  problem  asks  for  eopt{A,B), 
the  smallest  value  e  for  which  A  and  B  are  e-congruent.  It  is  possible  to 
use  any  algorithm  for  the  decision  problem  to  approximate  CoptiA,  B),  by  a 
search  procedure. 

Early  algorithms  for  the  e-congruence  decision  suffered  problem  from 
large  run-times.  A  concept  introduced  by  Stefan  Schirra  to  obtain  more 
efficient  algorithms  is  the  approximate  decision  algorithm.  A  decision  algo¬ 
rithm  is  {a, approximate,  if,  for  any  e  ^  \eopt{A,B)  —  a,eopi{A,B)  -)-  /3], 
it  correctly  answers  a  query,  and  for  e  €  [eopt[A,B)  —  a,eopt{A,  B)  -k  /5],  it 
either  answers  correctly  or  chooses  not  to  answer.  As  we  intuitively  believe 
that  it  is  more  difficult  to  test  for  e-congruence  if  e  is  close  to  £opt{A,  5),  we 
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would  be  willing  to  accept  an  algorithm  that  correctly  answers  queries  for 
values  of  e  not  near  £optiA,B),  but  sometimes  chooses  not  to  answer  when 
£  is  near  £opt{A,  B),  if  that  algorithm  provides  substantial  time  savings.  An 
(q,  /?)-approximate  algorithm  has  the  desirable  property  that  it  will  not  re¬ 
turn  an  incorrect  answer;  if  it  is  not  sure,  it  will  simply  say  that  it  does  not 
know  the  answer.  For  this  reason,  the  designation  “approximate”  is  perhaps 
a  misnomer;  an  alternate  title  for  an  algorithm  that  sometimes  chooses  not 
to  answer  is  “incomplete  decision  algorithm”  (as  opposed  to  a  “complete 
decision  algorithm,”  which  answers  all  queries). 

2  Generalized  approximate  algorithms 

While  “Approximate  decision  algorithms  for  point  set  congruence”  develops 
approximate  algorithms  for  the  original  formulation  of  e-congruencj,  this 
formulation  is  unnecessarily  constrictive,  which  motivates  the  desire  to  gen¬ 
eralize  results.  This  generalization  has  been  the  main  activity  of  the  project, 
and  the  principal  results  are  presented  in  the  paper,  “Generalized  approxi¬ 
mate  algorithms  for  point  set  congruence.”  While  this  paper  is  attached  to 
this  report,  we  will  summarize  briefly  its  content. 

First  of  all,  the  work  of  this  project  heis  dropped  an  assumption  of  the 
earlier  work:  the  point  sets  are  not  required  to  have  equal  cardinality.  This 
greatly  enlarges  its  applications.  Secondly,  the  work  examines  point  sets 
not  only  in  2-dimensional  space  but  in  3-dimensions  also;  it  even  begins  an 
examination  of  comparing  a  2-dimensional  set  to  a  3-dimensional  set. 

The  present  work  also  introduces  a  stronger  measure  for  congruency.  The 
original  decision  problem  merely  tells  whether  it  is  possible  to  match  the 
points  such  that  each  point  is  within  the  tolerance  distance;  the  present 
work  asks  for  the  maximum  number  of  pairs  that  can  be  placed  within  the 
tolerance  distance. 

This  notion  is  formalized  as  follows.  We  are  given  two  point  sets  A  and 
B  in  IR*^,  with  |A|  =  n,  |fi|  =  m,  n  <  m,  a  family  of  isometries  J,  and  a 
metric  dist{-,  •)  on  We  define  two  functions, 

^max  :R‘'’ -»  {2,  ...,n}  and  £opt  :  {2, . . .  ,n} -♦  R*', 

which  describe  the  congruence  between  A  and  B.  The  definitions  are  as 
follows. 
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•  We  say  that  A  and  B  are  {e,k)- congruent  if  there  exist  I  ^  I  and  an 
injection  (1-to-l  function)  £  :  A—*  B  such  that  dist{£{a),  /(a))  <  e  for 
at  least  k  points  a  G  >1. 

•  For  £  >  0,  kmaiis)  (or  kmaxi^.  A,  B))  is  the  maximum  value  k  such  that 
A  and  B  are  (e,/:)-congruent. 

•  For  k  €  {2, . . . ,  n},  Coptik)  (or  Coptik,  A,  B))  is  the  smallest  value  e  >  0 
such  that  A  and  B  are  (£,fc)-congruent. 

Combining  the  notions  behind  the  functions  k^ax  ^opt  leads  naturally  to  a 
structure  called  the  (£,A:)-map,  which  summarizes  the  congruency  relationship 
between  two  point  sets. 

The  investigator  has  generalized  the  definition  of  an  approximate  deci¬ 
sion  algorithm  to  correspond  to  the  generalized  definition  of  approximate 
congruence.  The  result  is  a  powerful  framework  for  studying  congruence  of 
point  sets.  The  investigator  has  developed  some  algorithms  for  the  various 
formulations  of  the  generalized  point  matching  question,  and  these  are  given 
in  the  attached  paper.  The  investigator  hopes  to  continue  to  improve  the 
algorithms  for  this  problem. 

3  Other  work 

While  the  work  of  this  project  has  been  primarily  concerned  with  point  set 
matching,  the  investigator  has  benefitted  from  the  support  of  the  grant  while 
pursuing  research  on  some  related  topics.  Point  matching  is  a  topic  in  com¬ 
putational  geometry,  since  it  asks  for  algorithmic  solutions  to  geometrical 
problems.  The  investigator  has  recently  worked  on  some  other  computa- 
tionaJ  geometry  topics  2is  well. 

The  investigator  hzis  studied  a  question  in  polygonal  visibility  known  as 
the  two-guard  problem.  Given  a  simple  polygon  P  with  vertices  s  and  <,  the 
straight  walk  problem  asks  whether  we  can  move  two  points  monotonically 
on  P  from  s  to  t,  one  clockwise  and  one  counterclockwise,  such  that  the 
points  are  always  co- visible.  In  the  counter  walk  problem,  both  points  move 
clockwise,  one  from  s  to  i  and  the  other  from  t  to  s.  Optimal  0(n)  con¬ 
structive  algorithms  for  both  problems  are  found.  The  results  are  obtained 
by  examining  the  structure  of  the  restrictions  placed  on  the  motion  of  the 
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two  points,  and  by  employing  properties  of  shortest  paths  and  shortest  path 
trees.  The  results  are  described  in  the  paper  “An  optimal  algorithm  for  the 
two-guard  problem,”  which  has  been  submitted  to  the  journal  International 
Journal  on  Computational  Geometry,  copies  of  the  paper  were  furnished  to 
the  A.R.O.  at  the  time  of  the  paper’s  submission  to  the  journal.  This  paper 
will  be  presented  at  the  9th  ACM  Symp.  on  Computational  Geometry,  May 
1993. 

The  investigator  has  researched  the  topic  of  t-spanners.  Let  K  be  a  set  of 
n  points  in  3-dimensional  Euclidean  space.  A  subgraph  of  the  complete  Eu¬ 
clidean  graph  is  a  f-spanner  if  for  any  u  and  v  in  V,  the  length  of  the  shortest 
path  from  u  to  t;  in  the  spanner  is  at  most  t  times  d{u,v).  The  investiga¬ 
tor,  with  Gautam  Das  and  Giri  Narasimhan,  has  shown  that  for  any  t  >  1, 
a  greedy  algorithm  produces  a  f-spanner  with  0{n)  edges,  and  total  edge 
weight  0(1)  •  wt{MST),  where  MST  is  a  minimum  spanning  tree  of  V.  This 
result  is  described  in  “Optimally  sparse  spanners  in  3-dimensional  Euclidean 
space,”  which  will  be  presented  at  the  9th  ACM  Symp.  on  Computational 
Geometry,  May  1993. 

The  investigator  recently  made  final  revisions  on  “Linear-time  algorithms 
for  weakly-monotone  polygons,”  which  has  been  accepted  in  revised  form  by 
the  journal  Computational  Geometry:  Theory  and  Applications.  The  paper 
introduces  a  new  class  of  simple  polygons  called  the  weakly-monotone  class. 
It  provides  a  lineax-time  algorithm  that  determines  whether  a  simple  polygon 
is  weakly-monotone,  and  a  brief,  linear-time  algorithm  that  triangulates  a 
weakly-monotone  polygon. 
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4  Papers 

“Generalized  approximate  algorithms  for  point  set  congruence,”  mcinuscript, 
1993. 

“Approximate  decision  algorithms  for  point  set  congruence,”  in  Proc.  of  the 
8th  ACM  Symposium  on  Computational  Geometry,  1992;  to  appear  in  Com- 
putational  Geometry:  Theory  and  Applications. 

“An  optimal  algorithm  for  the  two-guard  problem,”  to  be  presented  at  9th 
ACM  Symposium  on  Computational  Geometry,  San  Diego,  1993. 

“Optimally  sparse  spanners  in  3-dimensional  Euclidean  space,”  to  be  pre¬ 
sented  at  9th  ACM  Symposium  on  Computational  Geometry,  San  Diego, 
1993.  Co-Authors:  G.  Das  and  G.  Narasimhan. 

“Linear-time  algorithms  for  weaJdy-monotone  polygons,”  to  appear  in  Com¬ 
putational  Geometry:  Theory  and  Applications. 
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Generalized  Approximate  Algorithms 
for  Point  Set  Congruence 


Paul  J.  Heffernan  * 

Dept,  of  Mathematical  Sciences,  Memphis  State  University,  Memphis,  TN  38152 


Abstract 

We  address  the  question  of  determining  if  two  point  sets  are  congruent.  This 
task  consists  of  defining  the  concept  of  point  congruence,  and  developing  algorithms 
that  test  for  it.  We  introduce  the  (e,A:)-map,  a  device  which  pictoriaJly  represents  the 
degree  of  congruence  between  two  point  sets.  Point  set  congruence  has  been  studied  by 
previous  researchers,  but  we  feel  that  our  definitions  offer  a  more  general  and  powerful 
approach  to  the  problem.  By  using  the  paradigm  of  approximate  algorithms,  we  are 
able  to  construct  the  (£,fc)-map  efficiently. 


1  Introduction 

In  this  paper  we  address  the  following  question:  given  two  point  sets  in  IR*^,  what  do  we  mean 
when  we  say  that  they  are  “congruent,”  and  how  do  we  determine  congruency?  To  answer 
this  question,  we  must  define  congruency,  and  provide  efficient  algorithms  that  test  for  it. 
The  practical  motivation  behind  our  study  comes  from  computer  vision,  where  an  observed 
image  is  compared  to  a  hypothesized  model.  While  previous  researchers  have  formalized 
and  studied  this  problem,  we  feel  that  the  formulations  often  have  been  too  restrictive.  In 
particular,  an  elegant  formulation  has  been  introduced  by  [AMWW]  (also  [Ba])  and  studied 
further  by  [HS].  It  is  this  formulation,  which  we  call  e-congruence,  which  we  generalize  here, 
in  an  attempt  to  lessen  its  restrictions  while  preserving  its  strengths. 

Roughly  speaking,  the  e-congruence  formulation  is  as  follows.  The  input  is  assumed  to  be 
two  equal  cardinality,  planar  point  sets,  with  a  given  metric  (usually  L2  or  Too)  and  a  given 
family  of  isometries  (an  isometry  is  a  mapping  from  IR*^  to  itself  that  preserves  distances). 
There  are  two  problem  versions,  decision  and  optimization.  In  the  decision  problem,  one 
asks  for  a  given  value  £  >  0  whether  an  isometry  from  the  family  can  be  applied  to  one  point 
set  and  a  matching  found  between  the  sets  such  that  the  points  in  each  matched  pair  are 
within  distance  £  of  each  other.  In  the  optimization  problem,  one  asks  for  the  smallest  value 
£  that  returns  an  affirmative  answer  in  the  decision  problem. 

'Supported  in  part  by  the  U.S.  Army  Research  Office,  Grant  No.  DAAL03-92-G-0378 
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Formulation  of  the  model  as  a  point  set  is  appropriate  in  applications  where  the  mode! 
consists  of  a  number  of  small  pieces,  such  as  locations  on  a  map,  or  the  model  is  a  figure 
that  can  be  represented  by  its  key  boundary  points,  such  as  letters  of  an  alphabet.  The 
e-congruence  formulation  has  both  advantages  and  drawbacks.  One  of  its  main  strengths  is 
that  is  allows  for  the  real-world  possibility  of  noisy  input — that  is,  imprecisions  in  recording 
the  image  points — by  considering  a  pair  of  points  to  be  matched  if  they  are  merely  within  a 
tolerance  distance  e  of  each  other.  Every  point  of  both  the  model  and  image  is  important, 
since  each  point  must  be  matched  to  a  nearby  partner.  This  allows  us  to  distinguish  between 
models  that  vary  in  only  a  few  points.  On  the  other  hand,  if  even  one  image  point  encounters 
an  amount  of  noise  exceeding  the  tolerance  value,  we  may  falsely  answer  that  an  image 
does  not  resemble  a  model.  Another  concern  with  the  strict  definition  of  [AMWW]  is  the 
requirement  that  the  two  point  sets  have  equal  cardinality.  Among  the  imperfections  of 
image  registrations  are  missing  points,  which  fail  to  appear  on  the  image,  and  spurious 
points,  which  are  erroneously  introduced  into  the  image.  The  possibility  of  such  events 
renders  the  equal  cardinality  assumption  impractical.  Additionally,  in  [AMWW]  attention 
is  limited  to  2-dimensionai  point  sets. 

Our  goal  in  this  paper  is  to  abandon  the  rigidity  of  the  traditional  definition  of  e- 
congruency.  In  lieu  of  the  decision  problem,  we  propose  the  following  generalized  version; 
for  a  given  £  >  0,  find  the  largest  number  of  point  pairs  that  can  be  placed  within  distance 
£  of  each  other,  under  some  matching  and  isometry  from  the  family.  The  generalized  op¬ 
timization  problem  asks,  for  a  given  value  k,  to  find  the  smallest  e  for  which  the  decision 
problem  returns  a  value  of  at  least  k.  We  introduce  and  show  how  to  con.pute  a  structure 
called  the  {£,k)-map,  which  stores  the  complete  generalized  solution.  We  drop  the  assump¬ 
tion  that  the  point  sets  must  be  of  equal  cardinality,  and  we  expand  our  methods  to  higher 
dimensional  point  sets.  We  consider  a  variety  of  metrics  and  families  of  isometries.  We  even 
consider  the  case  of  projection  congruency,  in  which  a  2-dimensional  image  is  compared  to 
a  3-dimensional  model.  As  a  result  of  these  efforts,  we  build  a  more  powerful  framework  for 
congruency  of  point  sets. 

While  Alt,  et  al.  [AMWW]  give  algorithms  for  the  e-congruence  problem,  most  suffer 
from  high  worst-case  run-times.  A  way  to  obtain  efficient  algorithms  while  preserving  the 
formulation  is  through  the  use  of  approximate  algorithms  for  e-congruence,  first  introduced 
by  Schirra  [Scl,  Sc2].  Faster  approximate  algorithms  have  subsequently  been  developed  by 
fleffernan  and  Schirra  [HS].  In  order  to  develop  algorithms  for  our  generalized  definition  of 
point  set  congruence,  we  draw  on  the  techniques  of  [HS],  and  develop  a  body  of  approximate 
algorithms. 


2  Generalized  approximate  congruence 

In  the  point  set  congruence  problem,  we  are  given  sets  A  and  B.  We  wish  to  match  points 
of  A  to  points  of  B  so  that  as  many  pairs  as  possible  are  as  close  as  possible.  In  this  sense 
we  have  a  bicriteria  problem  with  a  trade-off  between  the  two  goals.  In  order  to  obtain 
workable  formulations,  we  can  fix  one  goal  while  optimizing  the  other.  Under  this  approach, 
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one  problem  consists  of  matching  the  points  in  order  to  maximize  the  number  of  close" 
pairs,  where  “close”  is  determined  by  a  fixed  tolerance  value.  The  corresponding  problem 
fixes  the  number  of  close  pairs  and  asks  for  the  minimum  possible  closeness  parameter.  In 
either  problem,  we  are  allowed  to  impose  an  isometry  from  a  given  family  on  one  of  the  point 
sets. 

We  formalize  these  notions  as  follows.  We  are  given  two  point  sets  A  and  B  in  (Rl^,  with 
|y4|  =  n,  )jB1  =  m,  n  <  m,  a  family  of  isometries  J,  and  a  metric  disi{-,  ■)  on  (R”^.  We  define 
two  functions, 

^’mar  :  IR'*’  — ►  {2, . . .  ,n}  and  :  {2,. . .  ,n}  (R'*', 

which  describe  the  congruence  between  A  and  B.  The  definitions  are  cis  follows. 

•  We  say  that  A  and  B  are  {e ,k)- congruent  if  there  exist  1  £  I  and  an  injection  (1-to-l 
function)  i  :  A  —*  B  such  that  dist{i{a)y  1(a))  <  e  for  at  least  k  points  a  £  A. 

•  For  e  >  0,  k^axi^)  ^max(£,  B))  is  the  maximum  value  k  such  that  A  and  B  are 

(£,A:)-congruent. 

•  For  k  £  {2,...,n},  eopt(k)  (or  eopt{k,A,B))  is  the  smallest  value  e  >  0  such  that  A 
and  B  are  (£,^)-congruent. 

Combining  the  notions  behind  the  functions  kmax  and  Copt  leads  naturally  to  a  structure  that 
we  call  the  (e,k)-map,  an  example  of  which  is  pictured  in  Figure  1.  The  (£,ilr)-map  plots 
k  against  £,  and  it  consists  of  both  a  YES  region  and  a  NO  region  of  ordered  pairs  (£,^‘)' 
such  that  A  and  B  are  (£,/:)-congruentif  and  only  if  {e,k)  is  in  the  YES  region.  The  {£,k)- 
map  allows  one  to  see  the  interaction  between  the  two  goals  of  obtaining  a  small  closeness 
parameter  and  a  large  matching.  One  wishes  for  a  point  in  the  upper-left  corner  of  the  map, 
while  the  YES  region  is  anchored  in  the  lower-right  corner.  The  boundary  between  the  two 
regions  is  a  step-wise  curve  from  lower-left  to  upper-right  (the  peculiarity  of  a  discontinuous 
curve  serving  a.s  a  boundary  occurs  because  k  assumes  discrete  values).  In  essense,  the 
(£,fc)-map  is  a  picture  of  the  congruence  between  A  and  B. 

The  (£,^)-map  offers  some  clear  advantages  over  the  traditional  definition  of  £-congruence. 
The  latter  deals  only  with  matchings  of  size  n,  and  therefore  limits  its  perspective  to  the  top 
line  of  the  (£,A:)-map.  Much  additional  information  is  contained  in  the  full  (£,/r)-map:  two 
pairs  of  point  sets  may  have  the  same  value  of  eopt{n)  but  a  sharp  difference  at  Soptin  —  i)  for 
a  small  integer  i.  Indeed,  large  recording  errors  for  a  few  image  points  can  increase  £op((n) 
greatly,  while  £opi(n  —  i),  where  i  equals  the  number  of  such  “outliers,”  gives  a  much  truer 
picture.  We  should  also  note  that  the  (e,k)-m&p  drops  the  assumption  of  equal  cardinality 
point  sets.  The  relaxation  allows  the  formulation  to  encompass  the  case  of  los  or  spurious 
points  in  the  image.  Once  such  points  exist,  the  idea  of  a  “perfect  matching”  between  the 
two  point  sets  becomes  less  meaningful,  and  the  concept  of  the  (£,A:)-map  more  appropriate. 
For  these  rea.sons  we  feel  that  the  (£,A:)-map  is  the  true  manner  in  which  to  think  about 
approximate  congruence. 
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3  Generalized  approximate  algorithms 

Approximate  decision  algorithms  were  introduced  in  [Scl,  Sc2],  and  expanded  in  [HS].  Here 
we  present  a  generalized  description  of  approximate  <dgorithms  suitable  for  (e,A:)*congruence. 
For  fixed  e  and  k,  the  decision  problem  asks  whether  input  point  sets  A  and  B  are 
congruent.  It  is  natural  to  believe  that  this  question  is  most  difficult  when  e  is  near  Soptik), 
so  we  will  allow  an  approximate  algorithm  not  to  answer  in  such  a  situation.  An 
approximate  decision  algorithm  for  (£,A:)-congruence  of  A  and  B  is  one  which  (1)  correctly 
determines  whether  A  and  B  are  (£,A:)-congruent  whenever  e  ^  [eop<(A;)  —  a,eopt[k)  +  /3], 
and  (2)  either  correctly  determines  (£,/!:)-congruence  or  chooses  not  to  answer  when  £  € 
[Soptik]  —  a,eopt{k)  +  0].  We  call  [£opi(fc)  —  oi,Sopt{k)  +  /?]  the  imprecision  interval.  A  major 
strength  of  approximate  decision  algorithms  is  that  they  never  return  an  incorrect  answer;  if 
a  particular  query  (£,k)  is  too  difficult,  the  algorithm  simply  says  that  it  is  unsure.  Previous 
approximate  decision  algorithms  have  been  defined  only  for  k  =  n. 

Earlier,  we  stated  two  optimization  problems:  maximizing  k  for  fixed  £  (i.e.  computing 
kmaxi^))  minimizing  e  for  fixed  k  (computing  €op<(fc))-  We  wish  to  extend  the  concept 
of  an  approximate  algorithm  to  these  optimization  problems. 

Let  us  first  consider  kmaxi^^)-  For  a  fixed  e,  the  ordered  pair  {e^k)  is  an  instance  of  the 
decision  problem,  for  k  =  2, ...  ,n.  For  all  k  <  kmax{^)i  A  and  B  are  (£,A:)-congruent,  while 
for  all  k  >  kmaxi^)  fhey  are  not.  Therefore,  a  (complete)  algorithm  for  the  A:,„oj;-optimization 
problem  partitions  the  set  {2, . . . ,  n}  of  values  for  k  into  two  sets:  {2, . . . ,  fcmai(e)}i  which 
we  call  the  YES  set,  and  {kmaxis)  +  1,  •  • . ,«},  the  NO  set  (one  of  these  sets  may  be  empty). 

An  approximate  algorithm  for  the  kmaxis)  problem  partitions  (2, . . . ,  n)  into  three  sets, 
{2, . . . ,  ki},  {fci  -f  1, . . . ,  fca},  and  {/:2  +  1, .  •  ■  ,n},  labelled  YES,  MAYBE,  and  NO,  respec¬ 
tively,  where  ki  <  k^axi^)  —  ^2  (one  or  two  of  these  sets  may  be  empty).  For  values  of  k  in 
the  MAYBE  set,  the  algorithm  is  unable  to  determine  whether  k  <  kmaxi^^)-  We  call  such 


4 


Figure  2:  The  (e,Ar,i/)-approximate  map 


an  algorithm  (a,  /?)-approximate  if,  for  every  k  in  the  MAYBE  set,  e  is  in  the  corresponding 
imprecision  interval,  i.e.  e  €  [£opt{k)  -  oi,£optik)  +  /?]. 

Now  let  us  consider  the  situation  of  a  fixed  value  k.  The  interval  IR'*’  of  possible  values 
for  £  is  divided  into  the  NO  subinterval  (0,eopt(A:))  and  the  YES  subinterval  [eop((^),oo); 
these  intervals  are  returned  by  a  (complete)  algorithm  for  tnv,  eopt{k)  problem  when  it  is 
given  input  k.  An  approximate  algorithm  for  the  £opt{k)  problem,  when  given  k,  returns 
intervals  (0,ex),  [£i,£2),  and  {£2,00),  labelled  NO,  MAYBE,  and  YES,  respectively,  where 
£1  <  £optik)  <  £2-  As  above,  the  MAYBE  set  represents  those  values  of  e  for  which  the 
algorithm  is  not  able  to  determine  (£,/;)-congruence.  Such  an  algorithm  is  (Q,^)-approximate 
if,  for  every  £  in  the  MAYBE  set,  e  is  in  the  corresponding  imprecision  interval;  that  is, 

£1,  £2  €  [£opt(^)  ~  £opt{k)  +  0\- 


4  Building  the  (£:,/c,i/)-approximate  map 

In  this  section  we  show  how  to  use  a  (7, 7)-approximate  algorithm  for  the  k^axi^)  problem 
in  order  to  build  a  {£,k,i/)- approximate  map  (see  Figure  2).  A  (£,fc,i/)-approximate  map 
is  similar  to  an  (£,fc)-map,  except  that  in  the  approximate  map,  £opt{k)  for  a  given  k  is 
represented,  not  as  a  single  point,  but  as  an  interval  of  length  no  more  than  u  that  contains 
£optik)-  The  (£,fc)-congruent  region  is  separated  from  the  non-congruent  region  not  by  a 
step-wise  curve,  but  by  a  region,  which  we  call  the  imprecision  region  (or,  the  MAYBE 
region). 

Our  approach  consists  of  a  two-part  search.  The  first  part,  which  we  call  Phase  I,  finds 
for  each  k  ^  {2, . . .  ,n}  an  interval  (c,2c)  that  contains  Coptik)-  In  Phase  II,  a  binary  search 
is  performed  on  each  interval  (c,  2c)  until  the  interval  size  is  at  most  i/. 
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Since  a  (7, 7)-approximate  algorithm  is  more  precise  for  small  values  of  7,  we  expect 
its  run-time  to  vary  inversely  with  7.  Actually,  we  will  see  that  the  run-times  for  our 
approximate  algorithms  for  the  kmaxi^)  problem  vary  directly  with  6/7,  so  we  attempt  to 
keep  this  value  small  when  constructing  our  (e,i)-map  algorithm.  We  will  see  that  the  time- 
complexity  of  Phase  II  dominates  that  of  Phase  I.  We  will  describe  the  general  (£,fc)-rnap 
estimation  method  now.  Later  we  will  give  (7,7)-approximate  algorithms  for  the  knaxi^} 
problem  in  various  dimensions  and  under  various  families  of  isometries  J,  and  after 
analyzing  each  such  algorithm  we  will  compute  the  time- complexity  of  the  corresponding 
(e,A:)-map  estimation. 

The  idea  behind  Phase  I  is  simple;  compute  tmar(£)>  initially  for  a  small  value  e,  then 
repeatedly  double  e  and  recompute  kmaxi^),  until  e  is  as  large  as  eopt{k).  At  each  step  we 
use  a  (7, 7)-approximate  algorithm,  where  we  set  7  =  e/2.  Let  us  first  describe  how  Phase  1 
works  for  a  single  value  k,  and  then  generalize  the  approach  for  all  k. 

We  begin  by  setting  e  =  i/,  and  then  alternating  the  steps  of  testing  for  (e,A;)-congruence 
by  using  a  (e/2,£/2)-approximate  algorithm,  and  doubling  e,  until  a  test  returns  an  answer 
of  YES  or  MAYBE.  (Note  that  an  answer  of  YES  (NO)  to  an  (e,A;)-congruence  test  implies 
that  e  >  €opt{k)  {e  <  Coptik)).)  At  this  point,  we  suspect  that  Coptik)  is  roughly  within  a 
factor  of  2  of  e,  since  the  test  for  (e/2,A:)-congruence  returned  a  NO  answer  while  the  one  for 
(e,A;)-congruence  returned  YES  or  MAYBE.  If  the  answer  for  (e, A:)- congruence  is  YES,  then 
we  realize  that  eopt{k)  €  (£/2,£),  which  is  an  interval  of  the  form  (c,  2c).  If  the  answer  is 
MAYBE,  we  see  that  Captik)  €  (e/2,3e/2),  since  we  used  a  (e/2,£/2)-approximate  algorithm 
to  test  for  (£,fc)-congruence.  This  interval  is  not  of  the  form  (c,2c),  but  we  can  trim  it  to 
this  form  through  one  additional  test  for  (e,A:)-congruence  with  7  =  e/4:  an  answer  of  YES, 
MAYBE,  or  NO,  respectively,  yields  an  interval  (e/2,  e),  (3e/4,5e/4),  or  (e,3£/2). 

There  is  one  special  case,  obtaining  a  YES  or  MAYBE  answer  on  the  initial  test,  when 
e  —  V.  Either  answer  implies  that  we  are  immediately  finished  not  only  with  Phase  I,  but 
with  Phase  II  as  well.  An  answer  of  YES  implies  eopt{k)  €  (0,  i/),  while  MAYBE  implies 
^opt{k)  €  (1//2, 3i//2);  thus,  in  either  case  we  have  found  an  interval  of  size  1/  that  contains 

^opt(^)' 

If  we  use  a  (7, 7)-approximate  algorithm  for  the  kmaxis)  problem,  then  we  can  perform 
Phase  I  simultaneously  for  all  values  k  =  1, . . .  ,n.  The  (7,7)-approximate  algorithm  returns 
an  answer  of  NO,  MAYBE,  or  YES  for  each  k,  so  we  alternate  performing  this  algorithm 
with  doubling  e  until  all  values  of  k  have  received  a  YES  or  MAYBE  answer  (fc  =  n  will  be 
the  la^t  to  receive  one  of  these  answers). 

Phase  II  is  quite  simply  a  binary  search.  Beginning  with  £opt(k)  €  (c,  2c),  we  proceed 
to  narrow  the  interval  to  width  i/.  A  step  begins  with  £opt(k)  €  (^,  u)-  We  test  for  (e,^)- 
congruence,  where  e  =  (^-j-u)/2,  and  7  =  (u  — £)/4.  With  a  YES  answer  we  set  u  *—  (^+u)/2 
and  repeat,  with  NO  we  set  €  (t ’  +  u)/2  and  repeat,  and  with  MAYBE,  our  choice  of  7 
allows  us  to  set  £  •*—  (3^  4-  u)/4,  u  *—  (£  +  3u)/4,  and  repeat. 
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5  Algorithms  for  kmaxi^) 

In  this  section  we  describe  (7,7)-approximate  algorithms  for  fcmor(£),  under  various  isom¬ 
etry  families.  We  will  use  the  Euclidean  metric,  although  our  algorithms  work  for  any  Lp 
metric.  Our  algorithms  are  based  largely  on  the  (7, 7)-approximate  decision  algorithms  for 
sT-congruence  of  [HS]. 

All  approximate  algorithms  employ  a  similar  strategy;  they  (1)  discretize  the  set  of 
isometries,  J,  into  a  finite  set  of  candidates,  C,  and  (2)  test  the  point  sets  for  congruence 
under  each  candidate.  The  methods  of  [HS]  employ  an  additional  tactic  of  perturbing  each 
point  in  A  and  B  in  order  to  impose  structure;  the  resulting  multisets  are  denoted  A*'  and 
B*.  This  method  consists  of  laying  a  lattice  onto  IR*^,  with  adjacent  lattice  points  distance 
A  apart,  where  A  is  proportional  to  7.  A*  and  are  obtained  by  moving  each  point  of  A 
and  B  to  the  nearest  lattice  point.  Note  that  no  point  is  perturbed  more  than  ■\/dA/2.  We 
will  see  later  how  working  with  A*  and  B*  instead  of  A  and  B  allows  an  improvement  in 
time-complexity. 

Each  algorithm  requires  that  the  following  conditions  hold  for  some  fixed  71,  72. 

Condition  1:  For  each  /  €  X  that  admits  (e,/:)-congTuence  for  some  A:  €  {2,. . .  ,n),  there 
exists  a  candidate  isometry  C  €  C  such  that  dist{C{a*)^I{a*))  <  71  for  all  a*  E  A*. 

Condition  2:  For  each  A:  <  n,  if  A*  and  B*  are  (e  —  72,  A:)- congruent  for  a  candidate 
isometry  C  €  C,  then  A  and  B  are  (e,A:)-congruent. 

Condition  3:  For  each  A:  <  n,  if  and  are  (e  -f  71  +  72,  A;)-congruent  for  no  candidate 
isometry  C  ^C,  then  A  and  B  are  not  (e,A;)-congruent. 

In  the  discussion  that  follows,  we  let  kmaxi/i,  A* ,  B*]  C)  represent  the  size  of  the  largest 
matching  between  A*  and  B*  with  tolerance  n  under  isometry  C.  Table  1  gives  an  approx¬ 
imate  algorithm  for  kmax{£iA,B),  assuming  we  have  constructed  C  and  have  a  procedure 
that  computes  k^,,{ti,A*,B*;C). 

[1]  Compute  set  of  candidates  C 

[2]  kyES  ^  1,  k^AYBE  *-  1 

[3]  forallCeCdo 

[4]  if  kmaxi^  +  7l  +  72,  A*,  B*\ C)  >  kMAYBE 

then  kMAYBE  kmaxis  +  7l  +  72,  B^\ C)  fi 

[5]  if  Ar„ax(e  -  72,  A*,B*;  C)  >  kyss 

then  kyES  fcmoi(e  -  72,  A*,  B*\ C)  fi  od; 

[6]  YES  ^  {2,...,  kyss} 

[7]  MAYBE  ♦— {AT'ss-f  l,...,fcAf4y be} 

[8]  NO  ^  {A:A/^yBE  + 1, •••,«} 

Table  1:  (7, 7)-approximate  algorithm  for  kmax{£,A,B) 
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Theorem  1  If  Conditions  2  and  3  are  satisfied,  then  the  algorithm  of  Table  1  is  a  (7,7)- 
approximate  algorithm  for  kmaxis,  A,  B),  where  7  =  71  +  272. 

In  our  algorithms,  72  =  2\/dX/2,  since  no  point  of  A  or  B  is  perturbed  more  than  vQA/2. 
We  will  also  choose  71  proportional  to  A  for  each  algorithm,  so  7  =  0(A). 

Implementation  of  an  algorithm  of  this  form  requires  (1)  a  method  for  constructing  C, 
and  (2)  a  method  for  computing  kmax{f^,A*,B*;C).  The  first  question  depends  on  the 
family  of  isometries  J,  and  will  be  explored  below.  The  second  question  leads  us  a  graph 
theoretic  formulation  and  solution. 

The  matching  problem  between  A*  and  B*  can  be  formulated  as  a  max-flow  problem. 
To  represent  graph  networks,  we  will  use  the  notation  {V,E,c),  where  V  and  E  are  the 
vertex  and  edge  sets,  respectively,  and  c  :  £  — »  IN  gives  the  edge  capacities.  Consider  the 
network 

G(C,  p.  A*,  B*)  =  ({5,  t}UUuV,E,  c) 

where  C  is  a  candidate  isometry,  U  =  represents  the  points  in  A*  and  V  — 

{uj, . . .  ,u„,}  the  points  in  B* , 

E  =  {(5,  w.)  I  1  <  i  <n}U  i(vj,t)  I  1  <J  <m}U  {{u.,u>)  |  disi(I(af),bf)  <  p}, 

and  c(e)  =  1  for  all  e  €  E.  The  max-flow  on  G  is  equal  to  k^axip,  A*,  B*\ C). 

At  this  point,  we  exploit  the  special  structure  of  the  perturbed  sets  and  B*,  in  order 

to  efficiently  compute  a  max-flow.  Feder  and  Motwani  [FM]  have  presented  a  technique  for 
improving  graph  algorithms  called  the  compression  graph,  which  is  well  suited  to  A*  and 
B*.  Let  =  {u,  I  a,  is  moved  onto  aj)  C  U  and  Vf  =  {vj  |  bj  is  moved  onto  6°}  C  V- 
If  C  moves  gridpoint  aj  the  /:i- neighborhood  of  b°,  this  generates  a  complete  bipartite 
subgraph  with  node  sets  and  Vf.  For  every  such  bipartite  clique,  we  remove  the  edges  in 

X  V°,  add  a  new  node  wm,  and  add  edges  1  €  C/^}  and  1  U;  € 

with  capacity  1  each.  Thereby  we  replace  If/^i  *  edges  by  jf/*  |  -f  [V)°|  new  edges.  W^e 

call  the  resulting  graph  B*). 

A  max-flow  in  G*^^p{C ,  p.  A* ,  B*)  corresponds  to  a  max-flow  in  G{C,  p,  A* ,  B*).  We 
can  show  that  Dinic’s  algorithm  computes  a  max-flow  on  p,  A*,  B*)  in  0{\/n  -f  m) 

phases,  each  in  time  proportional  to  the  new  number  of  edges.  The  reason  for  this  is  that 
Dinic’s  algorithm  takes  U(\/numfaer  of  nodes)  phases  in  a  simple  0-1  network.  A  node  is  called 
simple  if  it  has  indegree  1  or  outdegree  1,  and  a  network  is  called  simple  if  all  nodes  are  simple. 
Note  that  the  nodes  wu  are  not  simple  in  the  compression  of  a  graph;  however,  in  a  manner 
analogous  to  the  proof  on  the  number  of  phases  of  Dinic’s  algorithm  for  simple  0-1  networks 
in  [Me,  page  81],  it  can  be  shown  that  Dinic’s  algorithm  teJces  0(-v/numbcr  of  simple  nodes) 
phases  in  a  0-1  network,  if  no  edge  connects  non-simple  nodes. 

The  /i-neighborhood  of  any  point  contains  0((p/7)‘^)  grid  points  (since  7  =  0(A)).  Hence 
each  point  of  A*  U  B"^  contributes  to  0{{pf^Y)  bipartite  cliques.  Thus  the  number  of  edges 
in  graph  G*^^p{C,  p.  A* ,  B*)  is  0((n  -f  m){pl^Y),  which  yields  a  total  time  of  0((n  -|- 
my-^ipl'fY)  to  compute  kmax{fJi;A*,B*-,C),  The  graph  G*,^p  can  be  constructed  in  time 
0{(n  -f-  m)  log'^"*(n  -f  m))  using  standard  range-searching  techniques  [PSj. 
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Note  that  the  original  graph  G  is  a  simple  0*1  network,  for  which  a  max-flow  can  be 
computed  in  0{{n  +mY^^nm)  time.  Since  m  =  ©{«)  for  a  typical  problem  instance,  we  see 
that  the  use  of  the  compression  graph  that  is  enabled  by  the  structure  of  the  perturbed 
sets  A*  and  B*  reduces  the  run-time  in  terms  of  n. 

We  now  discuss  the  construction  of  the  set  of  candidates  C  for  various  isometry  groups. 
Throughout  this  section,  we  let  6  —  max{diam(y4),diam(B)}  be  the  larger  of  the  diameters 
of  the  sets  A  and  B.  For  each  isometry  group  J  we  must  find  a  set  of  candidates  C  that 
satisfies: 

Condition  1:  For  each  /  €  X  that  admits  (e,fc)-congruence  for  some  k  ^  {2, . . .  ,n},  there 
exists  a  candidate  isometry  C  such  that  dist{C{a*)^  I{a*))  <  71  for  all  a*  €  A* . 

If  J  equals  the  set  of  translations  on  then  let  Try  represent  the  translation  that  maps 
point  X  to  point  y.  By  fixing  point  x  and  choosing  as  candidates  all  translations  T^y  where 
2/  is  a  lattice  point,  we  satisfy  the  above  condition,  with  71  =  \/d\j2  {recall  that  each  point 
of  IR^'  is  no  more  than  distance  \/dXj2  from  the  nearest  lattice  point).  In  order  to  obtain 
a  finite  set  of  candidates,  we  note  that  if  a  candidate  C  admits  (£,t)-congruence  for  some 
k  >2,  then  C  maps  the  convex  hulls  of  A^  and  B^  within  a  distance  e  of  each  other.  Thus, 
a  collection  C  of  \C\  =  0{{{£  -f-  6)/Xy)  —  0({S/'fY)  candidates  suffices. 

Let  X  be  the  set  of  isometries  on  IR'^.  Since  an  isometry  is  an  affine  mapping,  it  can  be 
described  by  its  action  on  d  -i- 1  affinely  independent  points.  We  will  construct  a  collection  C 
that  satisfies  Condition  1,  where  each  candidate  isometry  is  a  set  of  d-f  1  affinely  independent 
points.  Fixing  the  first  point  is  equivalent  to  choosing  a  translation.  Once  the  first  point  is 
fixed,  each  remaining  point  is  free  to  move  on  a  d-dimensional  hypersphere  centered  at  the 
first  fixed  point,  and  after  d'  <  d  —  1  affinely  independent  points  have  been  fixed,  each  point 
not  in  the  affine  subspace  of  the  fixed  points  is  free  to  move  on  a  (d  —  d'  -f  1  )-dimensional 
hypersphere.  Once  d  affinely  independent  points  have  been  fixed,  each  point  not  in  the 
affine  subspace  can  assume  one  of  two  possible  values,  since  only  a  reflection  remains  to  be 
determined. 

We  construct  the  candidates  by  choosing  fixed  points  one  at  a  time.  For  the  first  point 
we  choose  0{{S/'<fY)  points  in  such  a  manner  that  every  translation  that  causes  the  convex 
hulls  of  A*  and  B*  to  intersect  has  its  fixed  point  within  distance  0(7)  of  the  fixed  point 
of  a  candidate.  When  d'  <  d  —  1  points  have  been  fixed,  the  set  of  potential  choices  for  the 
next  fixed  point  is  represented  by  all  points  on  a  hypersphere.  Since  no  free  point  moves 
on  a  hypersphere  of  radius  greater  than  S,  we  choose  0((^/7)^‘*"‘^  candidate  points 
on  the  hypersphere  of  radius  S  centered  at  the  first  fixed  point,  such  that  no  point  on 
this  hypersphere  is  more  than  distance  0(7)  from  a  candidate  point.  The  last  fixed  point 
represents  selecting  a  reflection  and  therefore  presents  only  two  choices.  If  the  dimension  d 
is  a  constant,  then  by  choosing  X  and  the  constant  in  the  0(7)  term  above  sufficiently  small, 
Condition  1  is  satisfied.  The  total  number  of  candidates  is 
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6  Time  complexity 

The  run-time  of  a  (7, 7)-approximate  algorithm  for  kmaz{£)  is  equal  to  0{T  ■  |(J|),  where 
\C\  equals  the  number  of  candidates,  and  T  the  time  to  test  each  candidate  of  C.  We  saw 
above  that  T  =  0{(n  + m)'  ®(e/7)‘^).  The  value  jC|  varies  with  the  isometry  family,  and  with 
the  above  observations  we  obtain  the  following  time-complexities  for  a  (7,7)-approximate 
algorithm  for  kmaxi^)-  0{{n  +  my-^{e/‘yY{6/'yY)  if  I  is  the  set  of  translations  in  (R*'  for 
A,  B  ^  (R*^,  and  0((n  -f  m)*  ®(e/7)‘'(i/7)‘*(‘^''’^l/^)  if  J  is  the  set  of  isometries  in  (R*^  for 

A,BeIR‘^. 

Let  us  examine  the  total  run-time  to  construct  a  (£,A:,i/)-approximate  map.  To  compute 
the  number  of  Phase  I  tests  for  a  given  k,  consider  that  we  test  for  (£,/:)-congruence  for  all 
values  £  =  2V,  e  =  0, . . . ,  j,  where£<,p<(A:)  €  This  implies;  =  ]og2i€optin)/i/) 

tests  until  we  receive  a  YES  or  MAYBE  answer  for  all  A:,  since  k  —  n  will  be  the  last  to 
receive  such  an  answer.  Also,  at  each  value  £  =  2V,  one  additional  test  may  have  to  be 
performed  because  certain  values  of  k  receive  for  the  first  time  there  an  answer  of  MAYBE, 
but  one  test  suffices  even  if  several  values  of  k  receive  a  MAYBE  at  the  same  value  £  =  2V. 
Thus,  the  total  number  of  values  of  £  tested  in  Phase  I  is  O(log(£opf(n)/t/));  since  7  =  £/2 
or  7  =  £/4  at  all  times  during  Phase  I,  the  total  time  is  0((n  +  log(£opt(n)/i/)), 

where  the  exponent  x  is  determined  by  the  family  of  isometries  J  amd  the  dimension  d. 

Phase  II,  for  a  given  k,  is  a  binary  search  on  an  interval  (c,  2c]  that  contains  £opt(k).  Each 
search  query  is  a  test  for  (£,A:)-congruence,  where  e  is  in  (c,  2c]  and  thus  approximately  equal 
to  £opt(k),  while  7  is  cut  in  half  after  each  query.  Therefore,  since  T  =  0((n  -I-  m)*  ®(£/7)‘^) 
and  \C\  varies  inversely  with  7  for  isometry  families  that  we  considered,  the  sum  of  all 
kmax{£)  tests  always  will  be  dominated  by  the  last  one,  for  which  7  is  approximately  the 
same  size  as  j/.  If  we  perform  Phase  II  over  all  n  veilues  of  fc,  we  get  a  total  time  of 
0{n{n  +  my-^{S/i'y{£opt{n)/i/Y),  with  x  as  above,  since  €opt{n)  >  £  for  every  test  value  £. 
We  see  that  Phase  II  dominates  Phase  I. 

In  summary,  implementing  Phase  I  and  Phase  II  together  builds  a  (£,fc,j/)-approximate 
map.  Using  the  (7, 7)-approximate  algorithms  for  kmaxi^)  that  we  have  described  yields  the 
following  time-bounds  for  constructing  (£,A:,i/)-approximate  maps: 

0{n(n  +  my-^{£opt(n)fi/Y{6/i/Y)  if  I  is  the  set  of  translations  in  IR'^,  for  A^B  €.  (R*^; 

0{n{n  +  rny  '^{€opt{n) / uy{S / if  J  is  all  isometries  in  IR*^,  for  A,B£  IR"^. 


7  Projective  congruence 

In  projective  congruence,  we  compare  a  2-dimensional  point  set  B  to  a  3-dimensional  set 
A.  Rather  than  choosing  just  an  isometry  and  a  matching  (between  the  elements  of  the 
point  sets),  finding  a  congruence  between  the  points  sets  now  consists  of  first  choosing  an 
orthogonal  projection  of  the  3-dimensional  set  A  toIR^,  and  then  choosing  an  isometry  on  the 
resulting  2-dimensional  set,  and  finally  a  matching.  We  develop  an  approximate  algorithm  for 
this  problem  by  approaching  it  in  a  similar  manner  as  before.  The  Phase  I-Phase  II  approach 
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motivates  a  version  of  the  kmax{£)  problem,  in  which  the  set  of  ail  possible  compositions  of 
projection  and  isometry  on  A  must  be  discretized  into  candidates  that  closely  approximate 
all  reasonable  choices.  Discretizing  the  candidates  resembles  the  same  task  for  isometries  in 
IR^,  except  that  the  projection  of  the  point  set  to  an  affine  space  of  one  lesser  degree  results  in 
the  loss  of  one  degree  of  freedom;  in  fixing  the  first  point  (which  determines  the  translation), 
we  realize  that  we  have  two  dimensions  of  freedom  and  not  three,  since  motion  parallel  to 
the  direction  of  projection  is  irrelevant.  The  total  number  of  candidates  on  the  last  test  of 
Phase  II  is  and  the  max-flow  problems  here  are  between  point  sets  in  IR^,  so  the 

total  run-time  is  0{n{n  -|-  m)^‘®(eop((n)/i/)^(6/t/)®). 


8  Conclusion 

We  have  generalized  the  concept  of  e-congruence  in  order  to  encompciss  a  greater  number  of 
situations;  our  generalized  definition  is  called  (e,A:)-congruence,  and  it  is  represented  pictori- 
ally  by  the  (£,/:)-map.  Because  we  desired  small  run-times,  we  chose  to  develop  approximate 
algorithms  that  test  for  (e,/:}-congruence,  using  many  of  the  techniques  of  [HS].  While 
we  did  not  discuss  complete  algorithms  for  (£,fc)-congruence,  such  algorithms  would  extend 
naturally  from  some  of  the  complete  algorithms  for  e-congruence  given  by  (AMWW). 

While  we  have  achieved  run-times  with  small  exponents  for  n  and  m,  it  is  troubling  that 
the  diameter  6  appears.  This  suggests  that  our  approximate  al  orithms  are  best  suited  for 
relatively  dense  point  sets.  Ironically,  for  e-congruence,  compRie  algorithms  perform  best 
on  point  sets  that  are  e-separated  [AMWW,  AKMSW]  (i.e.  no  two  points  of  the  same  set 
are  within  distance  e  of  each  other).  One  open  question  is  to  explore  which  types  of  point 
sets  are  best  suited  to  our  approximate  algorithms  and  which  to  complete  algorithms.  As  a 
final  open  question,  we  ask  if  there  are  ways  to  reduce  |C|  for  the  various  isometry  families 
we  studied,  since  any  reduction  leads  automatically  to  improved  algorithms  for  kmaxi^) 
the  (£,^,i/)-approximate  map. 
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